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3 Abstract. Neutrals sourced directly from Enceladus's plumes are initially 

4 confined to a dense neutral torus in Enceladus's orbit around Saturn. This 

5 neutral torus is redistributed by charge exchange, impact/photodissociation, 

6 and neutral-neutral collisions to produce Saturn's neutral clouds. Here we 

7 consider the former processes in greater detail than in previous studies. In 

8 the case of dissociation, models have assumed that OH is produced with a 

9 single speed of 1 km s~^, whereas laboratory measurements suggest a range 

10 of speeds between 1 and 1.6 km s~^. We show that the high-speed case in- 

11 creases dissociation's range of influence from 9 to 15 Rs. For charge exchange, 

12 we present a new modeling approach, where the ions are followed within a 

13 neutral background, whereas neutral cloud models are conventionally con- 

14 structed from the neutrals' point of view. This approach allows us to com- 

15 ment on the significance of the ions' gyrophase at the moment charge exchange 

16 occurs. Accounting for gyrophase: (1) has no consequence on the H2O cloud; 

17 (2) doubles the local density of OH at the orbit of Enceladus; and (3) de- 
is creases the oxygen densities at Enceladus's orbit by less than 10%. Finally, 

19 we consider velocity-dependent, as well as species-dependent cross sections 

20 and find that the oxygen cloud produced from charge exchange is spread out 

21 more than H2O, whereas the OH cloud is the most confined. 
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1. Introduction 

22 The Enceladus plumes directly produce a dense H2O torus centered on Enceladus's 

23 orbit, within which charge exchange and dissociation subsequently produce neutrals that 

24 either feed Saturn's extended neutral clouds, collide (absorb) with Saturn and its rings, 

25 or leave the system altogether on escape orbits. This paper is a report on the results 

26 of a sensitivity study of low- velocity charge exchange and dissociation within the neutral 

27 torus. 

28 Several decades before Cassini arrived at Saturn and the Enceladus water plumes were 

29 discovered [Hansen et al, 2006], neutral hydrogen was observed in Saturn's magneto- 

30 sphere, both from Earth [Weiser et al, 1977] and from Voyagers 1 and 2 [Shemansky and 

31 Hall, 1992]. Hydroxyl was later discovered by Shemansky et al. [1993] with HST, and 

32 more recently, Esposito et al. [2005] detected atomic oxygen. These observations coUec- 

33 tively hinted at the presence of a source of water, and models predicted its location to be 

34 near the orbit of Enceladus (c.f., Jurac et al. [2002]). 

35 After identifying the Enceladus plumes as the dominant source of the water-group 

36 neutrals (O, OH, H2O) — and indeed the co-rotating plasma itself via electron impact 

37 and photoionization [Young et al., 2005; Sittler et al., 2005, 2008] — researchers have been 

38 attempting to understand how neutrals are transported from Enceladus to 20 Saturn 

39 radii (Rs = 6 xlO^ cm) and beyond, as observed by Shemansky et al. [1993], Esposito 

40 et al. [2005], and most recently by Melin et al. [2009]. Early on, Jurac et al. [2002] 

41 mentioned the role of charge exchange in this inflation process. Johnson et al. [2006] later 

42 showed that if magnetospheric plasma is slowed sufficiently with respect to neutrals in 
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43 the Enceladus torus, then charge exchange produces a sufficient number of particles with 

44 velocities capable of spreading the dense H2O Enceladus torus into the cloud observed by 

45 Shemansky et al. [1993]. 

46 Farmer [2009] pointed out the importance of dipole-dipole interactions in collisions 

47 involving H2O molecules. She showed that collisions inside the dense Enceladus torus 

48 (parameterized by macroscopic viscosity) are alone sufficient for the creation of the ex- 

49 tended component of Saturn's neutral cloud. Cassidy and Johnson [2010] later argued 

50 that Farmer's fluid treatment is inappropriate for neutral-neutral collisions in the Ence- 

51 ladus torus, where the mean free path is on the order of the torus size itself. Instead, 

52 Cassidy and Johnson [2010] approached the problem with a direct simulation Monte Carlo 

53 (DSMC) model. Their model self-consistently included losses due to charge exchange, dis- 

54 sociation, and ionization, whereas Farmer [2009] accounted for losses to charge exchange 

55 and ionization by evolving the neutral cloud for the time scales (months to a few years) 

56 found in Sittler et al. [2008] . Both studies agree that neutral-neutral collisions are neces- 

57 sary for the inflation of Saturn's neutral cloud. 

58 Collisions between neutrals occur at a rate proportional to the square of the neutral 

59 density. Thus, where neutral densities peak near the orbit of Enceladus, neutral-neutral 

60 collisions occur more often than either charge exchange or dissociation, whereas near 6 

61 Rs, neutral densities drop and all three processes become comparable (see Fig. 3, Cassidy 

62 and Johnson [2010]). Models involving neutral collisions have recently been validated 

63 with Herschel observations by Hartogh et al. [2011], who attribute a warm and broadened 

64 Enceladus torus to heating via neutral-neutral collisions; the effect of these interactions 

65 should therefore be included in any attempt to fully model Saturn's neutral clouds. Nev- 
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66 ertheless, several first-order conclusions can be drawn by revisiting charge exchange and 

67 dissociation. 

68 Previous neutral cloud models approach charge exchange from the neutrals' point of 

69 view, whereas we follow the ion along its trajectory, thus allowing us to identify the 

70 gyrophase at which an ion undergoes charge exchange. We find that including the phase 

71 dependence doubles OH densities at the orbit of Enceladus, decreases oxygen density by 

72 < 10%, and has no effect on H2O (section 3.1). Also, the velocity-dependence of charge 

73 exchange varies by species. Previous studies {e.g., Johnson et al. [2006]; Cassidy and 

74 Johnson [2010]) have considered velocity-dependence, but have used a single cross section 

75 to represent all charge exchanges. We show in section 3.1 that symmetric reactions such as 

76 H2O -I- H2O+ — > H2O''" -I- H2O* (the asterisk identifies a neutral released with the speed of 

77 the reacting ion) tend to distribute neutrals closer to Saturn, while asymmetric exchanges 

78 such as H2O -I- 0+ — > H20''' -I- O* populate a more extended cloud, with less absorption 

79 on Saturn. 

80 With regard to dissociation, OH produced by impact/photodissociation of H2O has pre- 

81 viously been modeled with an initial speed of 1 km s~^ [Jurac and Richardson, 2005; Cas- 

82 sidy and Johnson, 2010], whereas recent laboratory measurements span 1 to 1.6 km s~^, 

83 depending on the molecule's internal energy [Wu and Chen, 1993; Makarov et al., 2004]. 

84 We model this parameter space and find that, relative to charge exchange, most OH found 

85 inside 9 Rs is produced by dissociation when OH is dissociated from H2O at 1 km s~^, 

86 with that location extended to 15 Rs when OH is dissociated from H2O at 1.6 km s~^ 

87 instead. 
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88 This paper is organized as follows. The model for the production of neutrals via dis- 

89 sociation and three illustrative charge exchanges is explained in section 2. Our results 

90 are found in section 3, followed by a discussion in section 4. The important points are 

91 summarized in section 5. 

2. Model 

92 We begin with a few words on nomenclature. The neutral torus in this paper pertains 

93 to the primary neutral torus (not plasma torus) supplied directly by Enceladus's plumes. 

94 The neutral clouds refer to the secondary neutrals produced from charge exchange and 

95 dissociation in the neutral torus. 

96 The production of Saturn's neutral cloud is modeled in two steps. We first construct 

97 a dense H2O torus from a plume positioned at Enceladus's south pole with specifications 

98 based on several Cassini Enceladus flybys {Smith et al. [2010]; see also Smith et al. [2004], 

99 Smith [2006], and Smith et al. [2007]). Secondary neutrals are then produced from the 

100 primary neutral torus by charge exchange and dissociation, some of which remain grav- 

101 itationally bound to Saturn and form the neutral clouds. On the basis that they spend 

102 most of the time outside the neutral torus and plasma sheet, we assume their lifetimes to 

103 be determined solely by photo-processes, though in section 3.3, we consider the effects of 

104 charge exchange and electron impact. 

2.1. Neutral torus model 

105 The neutral torus and Enceladus plume models are described in the following subsec- 

106 tions. 

107 2.1.1. Enceladus H2O torus 
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OS Our aim is to study the effects of several important reactions occurring in Enceladus's 
09 orbit. The primary neutral torus, fed directly by Enceladus, is produced in the model by 
evolving water molecules released from Enceladus into a dense neutral torus centered on 
Enceladus's orbit (3.95 Rg). The assumption is that all H2O is initially produced by a 
single plume at Enceladus's south pole. In reality, more than one plume has been observed 
[Porco et al., 2006], and researchers such as Saur et al. [2008] and Smith et al. [2010] have 

4 studied their signature on fly by observations. For our purposes, the detailed influence 

5 of multiple plumes can be neglected. The plume particles' radial speed distribution is 

6 prescribed as a one-dimensional Maxwellian with temperature T — 180 K [Spencer et al, 
2006; Hansen et al, 2006]: 



(V - Vbulk) 



2kT 

lis where Vbuik is the bulk speed, equal to 720 ms~^, 1.8x the thermal speed estimated by 

119 Smith et al. [2010] (I'therm = ^/'^kT/mn^o — 400 m s~^). Additionally, a raised cosine 

120 distribution is used to determine where the molecules are released: 



1 



1 + cos ( #-7r 



g(^0^ ^ ^ 9o ' V^o'VJ ^ < ^0 - 30° 
otherwise. 

121 Co-latitude theta is measured from Enceladus's south pole, and Oq — 30° is the plume half- 

122 width, based on INMS in situ observations [Smith et al, 2010]. We assume no azimuthal 

123 dependence. 

124 Enceladus's gravity is ignored since the escape velocity, v^sc = 240 m s~^, is greatly 

125 exceeded for most molecules (99%), where T;buik > Vtherm > Vesc', our results differ by less 

126 than one percent whether Enceladus's gravity is considered or not. Particles released from 

127 Enceladus are thus assumed to move on Keplerian orbits with respect to Saturn. Each 
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128 water molecule is allowed to orbit inside the torus for a period determined by the collective 

129 lifetimes against photodissociation, electron-impact dissociation, and charge exchange. To 

130 be clear, the molecules forming the neutral torus are subject to all of the losses stated, 

131 while the neutral cloud is subject to photo-processes only (section 2.2). Further reaction 

132 details are important for modeling plasma characteristics but only the timescales given 

133 below are required to model the neutral torus. 

134 Photodissociation: The photodissociation lifetime for H2O, Tphot — 9.1 x 10^ s, comes 

135 directly from Huehner and Carpenter [1979], scaled to Saturn's distance from the Sun. At 

136 peak solar activity, neutral abundances attributed to dissociation double [Jackman and 

137 Arndge, 2011]. 

138 Impact dissociation: In the case of impact dissociation, suprathermal (hot) electrons 

139 dominate [Fleshman et ai, 2010b]. Assuming the conditions near Enceladus's orbit apply 

140 throughout the neutral torus, we estimate the hot electron density and temperature to 

141 be 160 eV and 0.3 cm~^ from Fleshman et al. [2010b], which fit within a range of recent 

142 observations (cf.. Young et al. [2005], Sittler et al. [2008]). The rate coefficient for impact 

143 dissocation of water is found by convolving a{v)v with a 160 eV Maxwellian distribution 

144 to find Kimp = 1.5 X 10"^ cm^ s"^ {Fleshman et al. [2010b], table S9). Above 100 eV, 

145 Kimp is insensitive to temperature, making this estimate valid for a range of observations. 

146 Assuming the hot electron density, n^h, is constant over the neutral torus, the lifetime 

147 against impact dissociation is Timp ~ [^^imp^T-eh]"^ = 2.2 x 10^ s. Dissociation via thermal 

148 electrons — whose temperature and density are 2 eV and 60 cm~^ — is also expected, but 

149 such collisions occur 5x less often {Fleshman et al. [2010b], table S9) and are thus ignored 

150 here. 



DRAFT 



March 2, 2013, ll:02ain 



DRAFT 



FLESHMAN ET AL.: SATURN'S NEUTRAL CLOUDS X - 9 

151 Charge exchange: The following three reactions are included in this study: 

H2O + H2O+ ^ H2O+ + H2O* (3a) 
H2O + H2O+ ^ H3O+ + OH* (3b) 
H2O + O+ ^ H2O+ + O*. (3c) 

152 Other charge exchanges are important in the neutral torus, some of which involve ions 

153 reacting with secondary neutrals such as H, O, and OH. To model their effect properly, 

154 one would calculate these neutral densities as in a conventional time-dependent neutral 

155 cloud model. We estimate that including all such reactions would increase our estimates 

156 on neutral cloud densities by approximately a factor of two. The primary purpose of 

157 choosing this combination of reactions is to study three classes of charge exchanges, for 

158 which the collision frequency decreases with, increases with, or is independent of the 

159 relative speed of the reacting pair (reactions 3b, 3c, and 3a, respectively). We return to 

160 this point in section 2.2.2. 

161 An estimate of the charge exchange lifetime can be made by adding the rate co- 

162 efficients for reactions 3a-3c. Multiplying by the observed H2O+ and 0+ densities 

163 near the orbit of Enceladus (6 and 12 cm~^, Sittler et al. [2008]), we find Tchex — 
['«"xch^^H20+ + /«?xth^o+]"^ = [6.0 + 2.8]-^ X 10* s = 1.1 X 10^ s. The reaction rates 

165 are from Lishawa et al. [1990] and Albritton [1978] for /t^^h^ ^^"^ '^Sch; respectively. The 

166 hfetime of H2O against the sum of these processes is then 

-1 

= (1.1 + 4.5 + 0.91)-^ X 10^ s 

= 1.6 X lOS fti 20 days. (4) 



111 

+ + 



Tphot 
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167 These estimates can be compared to Fig. 3 of Cassidy and Johnson [2010]: for example, 

168 our lifetime against dissociation is 1.7 xlO^ s, while they use 7 x 10^ s near 4Rs. The 

169 discrepancy comes mostly from the impact dissociation timescale. Cassidy's dissociation 

170 rate was calculated using Schippers et al. [2008], wherein CAPS ELS data were fitted and 

171 extapolated down from 5.5 Rs, while our own estimate hinges on a hot electron density 

172 derived from our chemistry model [Fleshman et al., 2010b]. For charge exchange, we have 

173 a lifetime of 1.1 xlO^ s, and Cassidy and Johnson [2010] have a comparable 8 x 10^ s. 

174 Particles are created and tracked in each of our model runs, and the results are scaled 

175 to the number of water molecules in the real neutral torus. The total number is estimated 

176 from an assumed neutral source rate from Enceladus of M = 200 kg s^^ [Jurac and 

177 Richardson, 2005; Hansen et al., 2006, 2011] and lifetime, rtoms (Eq. 4): 

A^torus = Mrtorus/mHaO ~ 1.1 X 10^^ H2O molecules. (5) 

178 We bin and azimuthally average the results to find a 2-D density function, ntorus('", ^) 

179 (radius and latitude) , through which to introduce ions for charge exchange. This function 

180 also determines from where dissociated neutrals are produced. 

181 2.1.2. Plume model 

182 Before describing dissociation and charge exchange within the torus, we address a cal- 

183 culation with the purpose of comparing the neutral production near Enceladus with that 

184 from the entire torus. In doing so, we prescribe a plume whose density is consistent with 

185 Eqs. 1 and 2. In this case, the ambient neutral density can be ignored compared to the 

186 neutrals leaving the surface of Enceladus directly. 

187 The densities are determined everywhere by imposing integrated flux (J n(r, 9)v{r, 9)dA) 

188 and energy (mv^/2 — mM-EG/r) conservation at a given distance r. The picture can be 
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189 simplified, however, since most neutrals have at least twice Enceladus's escape velocity. 

190 The speeds are thus independent of r in the immediate vicinity of Enceladus. By equating 

191 the integrated flux at the surface of Enceladus to the same integral at another distance 

192 r > Re, we find a familiar dependence: 



193 The trailing exponential factor is imposed ad hoc to keep the total plume content finite, 

194 and reflects Saturn's influence as the molecules leave Enceladus. Consistent with Saur 

195 et al. [2008], is set at 4 times the Hill radius of 948 km. The angular dependence is 

196 consistent with our velocity distribution in the previous section (Eq. 1), 



197 which is normalized such that n(0) = no, where the plume strength uq is found to be 

198 5.9 X 10^ cm~^ by integrating n(^)vbuik over the area spanning the south pole of Enceladus 

199 from ^ = 0° to ^ = ^0 = 30°, and setting that result equal to the plume production rate 



201 In section 3, we compare the results of charge exchange with the plume (npiume) to that 

202 with the entire torus (ritorus)- 

2.2. Neutral cloud model 

203 Production of Saturn's neutral clouds entails following the neutrals produced by disso- 

204 elation and charge exchange occurring within the neutral torus. The treatment of each of 

205 these processes are described below. 

206 2.2.1. Dissociation 
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of M/mHaO = (200 kg s ^)/mH20 = 6.7 x 10^^ molecules per second. 
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207 The hydroxyl radical, OH, produced largely by dissociated H2O, has previously been 

208 modeled with a single speed of 1 km s~^ {i.e., Jurac and Richardson [2005], Cassidy and 

209 Johnson [2010]). Dissociated OH has been measured however with speeds between 1 and 

210 1.6 kms~^ [Wu and Chen, 1993; Makarov et al, 2004]. Here we bound this range by 

211 modeling the OH neutral clouds produced from an azimuthally-symmetric source (with 

212 respect to Saturn) with velocities drawn from Maxwellian distributions with tempera- 

213 tures T — ^moH'^mp' where the most probable speed, Vmp, is set to 1 and 1.6 km s~^, 

214 representing the low- and high-speed limits. 

215 The initial locations of the ejected OH are determined by the spatial distribution of 

216 neutrals in the Enceladus torus (ntomsl'^, 6*), section 2.1.1), and the directions of their 

217 release are chosen randomly and isotropically. The molecules orbit Saturn until they are 

218 photodissociated and removed from the system. 

219 By assuming a volume over which dissociations occur, the number of modeled OH 

220 molecules can be scaled to a realistic value. We take the volume to be a torus centered 

221 on Enceladus (3.95 Rs), with a minor radius of 1 Rs: 

y f« 27r(4 Rs) (2 Rs)^ = 2 X 10^^ cm^ (8) 

222 For impact dissociation, we then expect a contribution of 

^doud = ^impTp^hotV^ = 2.8 X 10^^ OH molecules, (9) 

223 where kuap = 7.9x 10~^ cm~^ s~^ is the rate (per volume) of impact dissociations occurring 

224 between suprathermal electrons and H2O molecules in the torus [Fleshman et al., 2010b] 

225 and r^Q^. = 1.8 x 10^ s is the photodissociation lifetime of OH at Saturn [Huebner and 

226 Carpenter, 1979]. The number of OH molecules produced by photodissociation in the 
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227 torus is similarly given by 

^Sud = W^p°hot^ = 7.6 X 10=^3 OH molecules, (10) 

228 where /Cphot — 2.1 x 10~^ cm~^ s~^ is the rate (per volume) of H2O photodissociations 

229 occurring in the Enceladus torus {Fleshman et al. [2010b], Table S9). The total abundance 

230 attributed to dissociation is then given by the sum of Eqs. 9 and 10. Cassidy and Johnson 

231 [2010] constrained their study with HST observations [Melin et al, 2009] and found a 

232 similar OH content (see comparison in Fig. 10c, this paper). 

233 That neutral production by photo- and impact dissociation are comparable in magni- 

234 tude is itself noteworthy. This condition is not shared by systems with hotter and denser 

235 plasma. For example, electron impact dissociation and ionization dominate over photon- 

236 driven processes in Jupiter's lo torus, where the plasma is warmer where the pick-up 

237 energies are four times higher than at Enceladus [Delamere et al., 2007; Fleshman et al., 

238 2010b]. We also note that unhke with lo, long neutral hfetimes in the Enceladus neu- 

239 tral torus inhibit the response of Saturn's neutral clouds to short-term plume variability, 

240 though variabihty on the order of months has been studied by Smith et al. [2010]. 

241 2.2.2. Charge exchange 

242 We now describe the model for producing and following neutrals from charge exchange. 

243 Cassidy and Johnson [2010] and Jurac and Richardson [2005] also considered velocity- 

244 dependent charge exchange, but unlike these previous studies, we capture the gyrophase 

245 at which the reactions occur by following ions along their trajectories (section 3.1). We 

246 also prescribe cross sections specific to each reaction, being particularly interested in the 

247 effects of low- velocity charge exchange. 
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248 At very high speeds, the cross sections go to zero for all charge exchanges [Johnson, 

249 1990]. At low relative velocities, however (few km s"^), the details of the coUision are 

250 determined by the nature of the reacting species. If the reactants and products are 

251 identical, apart from an electron (i.e., H2O + H2O+ — > H2O+ + H2O*), the reaction is 

252 termed resonant, or symmetric, and the cross sections grow as the inverse of the relative 

253 speed. If the reactants differ, as with 1120 + 0"^ — > H2O+ + O*, the cross sections are 

254 likely to vanish at low speeds [Rapp and Francis, 1962] — the difference being that the 

255 energy of the electron configurations is unchanged for symmetric-type charge exchanges 

256 [Johnson, 1990]. Neutrals produced from resonant charge exchange therefore tend to 

257 have lower velocities than do neutrals produced from non-resonant (asymmetric) charge 

258 exchange. This is a key point central to much of our discussion in section 3. 

259 Individual ions are followed as they traverse the neutral torus (section 2.1). This ap- 

260 proach allows their gyrophase to be determined the instant that charge exchanges occur 

261 (see Fig. 1). The implicit assumption is that the collision is elastic, and that the neutral 

262 product has an initial velocity given by the ion velocity just before the exchange takes 

263 place. 

264 The ions are introduced into the model from two Maxwellian speed distributions. 




.2 1 





(speeds parallel to B) 




265 



with a temperature anisotropy of 



kTi 27 eV 



(13) 



A;7|| 5.4 eV 



5 
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266 for both O"*" and H2O"'" [Sittler et ai, 2008]. The perpendicular temperature is derived 

267 from the pick-up ion vefocity at the orbit of Enceladus, determined from CAPS data by 

268 Wilson et al. [2009] [kT±_ = \m-si^+{y^ — fRep)^]- The ions also rotate around a guiding 

269 center (field line) moving at — 18 km s~^ in a frame rotating with the neutrals. 

270 For the component of our study aimed at estimating local neutral production (section 

271 2.1.2), ions passing near Enceladus are diverted (treating Enceladus as a rigid cyhnder) 

272 and are slowed to 10% of the ambient flow speed to account for the effects of mass-loading 

273 (see Fleshman et al. [2010a]). 

274 Time steps are taken at less than 1% of an ion's gyroperiod: 

At^Rx Tgyro ^Rx (14) 

275 where R is a random number between and 0.01, Tgyro is the ion's gyroperiod (3.6 s for 

276 H2O"'"), B = 325 nT, and q and m are the charge and mass, respectively, of the reacting ion. 

277 Such resolution is necessary in order to capture the significance of the energy dependence 

278 at low relative speeds. After each time step, the collision frequency u is calculated from 

i/(r, 9, t^rel) = n{r, 9)a{Vrel)Vreh (15) 

279 where n{r,9) is the local H2O density (section 2.1), v-^ei is the relative velocity between 

280 the reacting ion and neutral, and cr(t'rei) is the velocity-dependent cross section. Poisson 

281 statistics are used to test the likelihood of one or more reactions having occurred within 

282 At. If exp(— i/At) is less than a second random number between and 1, then a reaction 

283 occurs. The possibility of multiple reactions occurring over At is taken into account, but 

284 it is neglectable (appendix A). 
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285 As with OH produced by dissociation (section 2.2.1), neutrals produced by charge ex- 

286 change are followed under the influence of Saturn's gravity until they are photodissociated 

287 or photoionized. Their initial location and velocity are taken to be that of the reacting 

288 ion, pre-transfer. 

289 The model runs are centered on the orbit of Enceladus spanning 10 Re in the direction of 

290 corotation (Re = 250 km — radius of Enceladus) and ± 120 Re (0.5 Rs) in both the radial 

291 and z directions to adequately sample the H2O torus (section 2.1.1). Ions are introduced 

292 into the model on the upstream boundary, and their guiding centers flow downstream at 

293 a speed I'piasma = 18 km s~^ relative to the neutrals. Their starting location in (r, z) is 

294 chosen randomly. 

295 Scaling: The neutral clouds formed via charge exchange are done so in our model by 

296 following a relatively small number of ions, and must thus be scaled to facilitate compar- 

297 ison with observations and other models. The number of neutrals in our modeled clouds 

298 have been scaled by accounting for the following. First, the number of representative ions 

299 used to produce the neutral clouds via charge exchange falls short of, and must be scaled 

300 to, the number of ions present in the actual plasma torus, n\o^ ■ The volume of the 

301 plasma torus, V ^ is given in Eq. 8, and n = 12 and 6 cm~^ for O"*" and H2O+, respectively 

302 [sutler et al., 2008]. Second, we have argued that photo-processes are more likely to occur 

303 than either charge exchange or electron-impact processes throughout the neutral clouds 

304 with the exception of very near the neutral torus. In keeping with this assumption, the 

305 plasma torus thus feeds the extended neutral clouds via charge exchange for a photodisso- 

306 ciation (photoionization in the case of oxygen) time scale before equilibrium of the neutral 

307 cloud is achieved: Tphot = 14, 0.6, 0.3 years for O, OH, and H2O, respectively [Fleshman 
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308 et ai, 2010b]. Our model runs followed 10^ ions for 100 seconds, and the resulting neutral 

309 clouds were scaled as described. 

3. Results 

310 In the following sections, we present and discuss the neutral clouds resulting from dis- 

311 sociation and charge exchange in our model. 

3.1. Chcirge exchange 

312 In the neutrals' reference frame, ions oscillate between km s~^ and twice the local 

313 pick-up speed {v^ ~ 18 km s~^) due to gyro-motion. A cartoon of this can be seen in Fig. 

314 1, where v^ei ~ at the cusp of the ion trajectory and reaches a maximum of v^ei ~ 

315 along the flow direction. Shown are several trajectories for which v± is either less than, 

316 greater than, or approximately equal to the bulk flow velocity. The neutrals formed via 

317 charge exchange follow the trajectories indicated in red. 

318 The velocity dependence of reactions 3a-3c are determined by the details of the reacting 

319 species [Johnson, 1990]. Essentially, symmetric charge exchanges have cross sections that 

320 increase monotonically with decreasing velocity, whereas cross sections for asymmetric 

321 exchange peak and then vanish at low relative speeds. The implication is that symmetric 

322 exchanges produce lower velocity neutrals and a more compact neutral cloud than do 

323 asymmetric reactions. 

324 With symmetric charge exchange, the cross sections go as Vj-eC^, so that the collision 

325 frequency (nav) is independent of f , as with reaction 3a, whereas asymmetric exchanges 

326 are defined by cross sections (and collision frequencies) which tend rapidly toward zero at 

327 low relative velocities (~ v^ei^, Rapp and Francis [1962]). 
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328 The cross sections (10~^^ cm^) used in this paper to study reactions 3a-3c plotted in 

329 Fig. 2a are given by 

Reaction 3a : H2O + H2O+ ^ H2O+ + H2O* 

(7H,o = 38£;-°-5 (16a) 
Reaction 3b : H2O + H2O+ ^ H3O+ + OH* 

aon = 38£;-r - 0.39exp ( %^Vl (16b) 
Reaction 3c : H2O + 0+ ^ H2O+ + O* 
ao = 69£;-°-29 + 30exp 



1 /E,,,-65^ ' 



2V 18 ' ■ ^'''^ 

330 The Gaussian terms in Eqs. 16b and 16c account for downward and upward trends in the 

331 associated data sets near 30 km s~^, but have little consequence on the neutral cloud, 

332 given that most bound particles are produced at lower velocities. 

333 Symmetric exchanges occur between like species by definition, although unlike species 

334 also exhibit symmetric behavior on occasion. Therefore, we explore several hypothetical 

335 behaviors for the OH*-producing reaction 3b at low energies. This test is separate from, 

336 but related to, the comparison between reactions 3a-3c themselves, and it motivates the 

337 point that both high and low energy behaviors have an important effect on the neutral 

338 cloud. With o-^^^''ap°i3.ted^ hscve extropolated the best-fit curve (Eq. 16b) to the lowest 

339 energies. Symmetric and asymmetric behaviors are explored with ^nd (j^^y™™^*'"''^ 

340 [Rapp and Francis, 1962; Johnson, 1990]. jg the same as except that 

341 below 1.5 eV, a^y^""^^"" = SOE'^-^ x 10"^^ cm^. Notice that a similar energy dependence 

342 also applies to Eq. 16a, consistent with symmetric charge exchange. Q-^y^^^^"'^ jg the same 

343 as a^^^""'^'^^ except that below 1.5 eV, cr^y--^*"'^ = UEl, x lO'^^ ^^^^2 Although it could 
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344 be argued that (jg'™™^*''^'^ better fits the data if the two measurements at 2 eV are ignored, 

345 our results for reaction 3b were obtained with (j^^^p°^^^^'^ unless noted otherwise. We will 

346 discuss the implications of choosing q-^^^p"^^^'^'^ over (jg'™™''*''"^ ^nd (j^sy™™^*'''^ shortly. 

347 The collision frequencies (nav) are plotted in Fig. 2b for a given neutral density — in 

348 this case for riHjO — 10^ cm~^. The collision frequency for oxygen increases with relative 

349 speed, while it is constant for water, and peaks at low velocities for OH. The significance 

350 is that the oxygen cloud tends to be more extended than either the OH or H2O clouds. 

351 The average collision frequency is also much higher for oxygen (x 10) than for either OH 

352 or H2O, resulting in greater oxygen abundance. 

353 The equatorial neutral cloud densities resulting from reactions 3a-3c are plotted in 

354 Fig. 3. Only neutrals produced from charge exchange are shown; neither the Enceladus 

355 neutral torus, nor the neutrals produced via dissociation have been included. Oxygen is 

356 two orders of magnitude more abundant than either OH or H2O because of the higher rate 

357 of production, but also because oxygen has a longer lifetime against photoionization than 

358 either OH or H2O have against photodissociation. Unlike Cassidy and Johnson [2010], 

359 dissociated neutrals from the latter processes are not tracked in our model. Beyond the 

360 scope of the present study, this additional heating source would serve to further inflate the 

361 oxygen and OH clouds. Fig. 3b is the same as 3a, except that the profiles are normalized 

362 to the peak density at the orbit of Enceladus. The oxygen cloud is seen to be the most 

363 extended, followed by water, and finally by OH, with an order of magnitude separating 

364 the three species at 20 Rs. 

365 The effects of low-velocity charge exchange are shown in Fig. 4. In Fig. 4a, we see that 

366 the peak density (as well as the total neutral cloud content) is the highest with 
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I 1 I -J- J- 1 J J j-u -j-U -j-U symmetric asymmetric 

367 because more low-velocity neutrals are produced than with either cTqjj or cToh 

368 Conversely, fewer low-velocity neutrals are available to populate the region near Ence- 

II ) 1 •! -ii asymmetric i j i -.i extrapolated symmetric cj- -t J 

369 ladus s orbit with ctq^ when compared to either ctqjj ^ or ctqjj . btated 

370 another way, yields a neutral cloud with the steepest slope, and ^Jig 

371 shallowest. Fig. 4b is identical to Fig. 4a, apart from normalization. In this case, the slope 

372 of the density profile should not be confused with the effect of inflating (spreading) the 

373 OH cloud. It should be viewed, rather, as the enhancement or depletion of low velocity 

374 neutrals to fill the region inside of 10 Rs. In other words, neutrals beyond 10 Rs are 

375 mostly formed in charge exchanges at high velocities, for which all (Joh converge to the 

376 same curve (Fig. 2). 

377 We have assumed to this point that the plasma is sub-corotating in Enceladus's orbit 

378 (18 km s~^, Wilson et al. [2009]). One might expect, however, that the neutral clould 

379 would be affected in a measurable way if instead, the plasma corotates at 26 km s~^. 

380 The H2O cloud would be least affected, given that the collision frequency of reaction 3a is 

381 independent of speed (Fig. 2b), but what about reactions such as 3b and 3c, whose collision 

382 frequencies are velocity-dependent? Increasing the plasma speed amounts to shifting the 

383 spread of ion velocities in Fig. 2 to the right, which would on average increase the speed 

384 of the neutral products. This is indeed the case, and in such a test where we increased the 

385 plasma speed from 18 to 26 km s^^, the oxygen cloud increased in abundance and became 

386 even more extended. The OH cloud also expanded somewhat, but decreased in total 

387 abundance. Unfortunately, the differences were less than 10% in both the slope of the 

388 distribution and in total oxygen abundance, suggesting that neutral cloud observations 

389 are in this way unlikely to predict plasma speeds in the torus. 
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390 3.1.1. Neutral cloud sources: plume vs. neutral torus 

391 We described in section 2.1 the production in our model of the neutral H2O torus from 

392 the Enceladus plumes. The plumes themselves have also been prescribed as a separate 

393 background density (ripiume, section 2.1.2) so that we can compare charge exchange occur- 

394 ring throughout the neutral torus to that occurring only within the Enceladus plumes. 

395 The results are shown in Fig. 5, where we have plotted the oxygen clouds produced 

396 from charge exchange within both the Enceladus plumes (local) and the entire neutral 

397 torus (global). The results are for reaction 3c, but the same test with reactions 3a and 

398 3b produces similar results. Immediately noticable is that the local production is fa 0.1% 

399 of the overall neutral production. The torus's dominance of neutral production can be 

400 explained as follows. First, the volume of the torus where reactions are occurring can 

401 be estimated as 27r(4Rs)(0.2Rs)^, where O.lRs is roughly the torus's scale height. The 

402 volume of the plume can be estimated from Eq. 6, where the dimensions are on the order 

403 of a cylinder with width 2Re and height Hr ~ 16Re. Dividing these volumes gives roughly 

404 250(Re/Rs)^ ~ 10~^. Further, the collision frequencies are proportional to the neutral 

405 density, which in the plume are on the order of 10'' cm~^, whereas typical torus densities 

406 are 10^ cm~^, making collisions in the plume 100 x more frequent per volume than in the 

407 torus. All told, the ratio of the volumes (10~^) combined with the ratio of densities (10^) 

408 explain the local-to-global neutral production ratio of 10~^ shown in Fig. 5a. A similar 

409 pattern has been shown to exist at Jupiter by Bagenal [1997] and Dols et al. [2008] , where 

410 the majority of plasma is produced throughout the neutral torus, rather than near the 

411 interaction at lo itself. 
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412 The slopes of the neutral clouds from the plume and torus are most easily compared 

413 in Fig. 5b, in which the density profiles have been normalized. The local source produces 

414 a more confined neutral cloud because the ions from which they originate have been 

415 slowed near the plume to account for the effect of mass-loading (Fleshman et al. [2010a]). 

416 Nevertheless, such a signature would be difficult to untangle in the data since global 

417 exceeds local production so overwhelmingly. 

3.2. Dissociation 

418 A major component of the OH cloud is produced by dissociation within the neutral 

419 torus, whereby the initial velocities of the OH products range from 1 to 1.6 km s^^ [Wu 

420 and Chen, 1993; Makarov et al, 2004]. In Fig. 6, the clouds resulting from the high- and 

421 low-speed cases are plotted along with the result from velocity-dependent charge exchange 

422 in section 3.1. First note that dissociation contributes lOOx more OH than does charge 

423 exchange at the orbit of Enceladus (4 Rs); the total cloud mass is almost 100 x greater 

424 as well. Second, dissociation dominates over charge exchange from the Enceladus torus 

425 out to 9 and 15 Rs in the low- and high-speed cases, respectively. The OH cloud content 

426 will only be marginally affected by variable solar activity [Jackman and Arridge, 2011], 

427 given that impact dissociation contributes 4x more neutrals than does photodissociation, 

428 by virtue of the respective reaction rates (section 2.2.1). In both cases, few neutrals are 

429 absorbed by the rings, and even less by Saturn itself. The same is not true of charge 

430 exchange, where ~ 50% of the neutrals are absorbed by Saturn (section 3.3). 

431 Fig. 7c is a two-dimensional version of Fig. 6, where the dissociation results have been 

432 averaged and added to the results from charge exchange. Saturn is at the left, and the 

433 Enceladus 's orbit is located on the equator at 4 Rs. In addition to being confined radially. 



DRAFT 



March 2, 2013, ll:02ain 



DRAFT 



FLESHMAN ET AL.: SATURN'S NEUTRAL CLOUDS X - 23 

434 the dissociated neutrals are also bound tightly to the equator, while neutrals from charge 

435 exchange tenuously fill the magnetosphere. 

436 Fig. 7a shows the hydrogen cloud that accompanies the dissociated OH clouds 

437 (H20 + e, 7 — > OH* + H*). To conserve momentum, the hydrogen atoms have 17x the 

438 speed of the dissociated OH molecules, and thus range between 17 and 27 km s~^, with 

439 a relatively large, diffuse neutral cloud. Shown is the result for the low-speed case, which 

440 produces more bound particles and thus a more substantial neutral cloud. Charge ex- 

441 change from reactions such as H2O + H+ — > H2O"'" + H* are also responsible for H-cloud 

442 production, and deserve attention in future studies. 

3.3. Fates of neutral atoms and molecules 

443 In our model, neutrals created by dissociation and charge exchange are eventually either 

444 absorbed by Saturn, escape the system, or orbit until they are destroyed (ionized) by 

445 photons. In Fig. 8a the fates for each species are given by percentage. In the case of 

446 hydrogen, the results are from the dissociation model, described in section 3.2. The 

447 enormous amount of escape (84%) is due to the high velocities (~ 17 km s~^) with which 

448 hydrogen is created following H2O dissociation, and the 8% absorption is largely comprised 

449 of hydrogen which would escape the system otherwise. 

450 Oxygen is produced purely from charge exchange in our model (reaction 3c). About one- 

451 half escapes, one-third is absorbed, and the remaining 13% contributes to the neutral cloud 

452 before being photoionized. Water is also produced purely by charge exchange (reaction 

453 3a) with 18% contributing to the neutral cloud. Percentage-wise, more water is absorbed 

454 than oxygen because oxygen is produced with higher speeds and generally larger orbits 

455 (section 3.1). 
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456 The fate of OH is dominated by dissociation: 96% feed the neutral cloud (ultimately 

457 ionized), 4% are absorbed, and virtually none escape. The reason for the large percentage 

458 of bound and unabsorbed neutrals is that dissociated OH has a velocity spread of 1 to 

459 1.6 km s"^ in the neutral frame, compared to the escape speed of f=:i 5 km s"^ in the same 

460 frame. Looking only at OH produced by charge exchange (minor compared to dissocia- 

461 tion), 58% are absorbed, 23% supply the neutral cloud, and 20% escape. Compared to 

462 H2O, an even greater percentage of charge-exchanged OH is absorbed because the cross 

463 sections favor production of low- velocity OH molecules (Fig. 2b). 

464 The production of oxygen via dissociation of H2O has been ignored in this paper on 

465 the grounds that, unlike OH, oxygen is largely produced by charge exchange. The cross 

466 section for oxygen-producing charge exchange is an order of magnitude higher than that 

467 for the OH-producing reaction near the plasma flow speed of Vpiasma = 18 km s~^ (Fig- 

468 2a), while the photodissociation rates are an order of magnitude smaller [Huebner and 

469 Carpenter, 1979]. We estimate that including oxygen produced from dissociation would 

470 increase the total oxygen cloud content by less than 20%. 

471 Charge exchange and dissociation play a large role in creating Saturn's neutral clouds 

472 from the plume-fed neutral torus. The reactions we have included have been chosen to 

473 demonstrate the effects of low velocity charge exchange and dissociation, but they are 

474 also among the most important. The neutral cloud densities presented in this paper are 

475 expected to undershoot the results from models which include the additional reactions 

476 found in Fig. 3, of Fleshman et al. [2010b] by no more than a factor of two. With this 

477 caveat in mind, we now compare the present results with several other recent models. 

478 3.3.1. Comparison with other models 
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479 Fig. 8b: JOG is the work of Johnson et al. [2006], where they also investigated the neutral 

480 clouds created from low-velocity charge exchange in the stagnated flows in Enceladus's 

481 orbit. Fig. 8b: J07 is from Jurac and Richardson [2007], where the authors were primarily 

482 interested in the interaction between the neutral cloud and Saturn's rings. The most 

483 recent model comes from Cassidy and Johnson [2010] (CIO), where they investigated the 

484 spreading of the neutral cloud from neutral-neutral collisions. 

485 To compare with these studies, we first had to weight our H, O, OH, and H2O clouds. 

486 We did this for two limiting cases. In the first case (rphot. Fig. 8b), we assume, as we 

487 have thus far, that the neutral clouds evolve until destroyed by either photoionization or 

488 photodissociation: H, O, OH, H2O = 40, 14, 0.6, 0.3 years, respectively. These lifetimes 

489 yield an upper limit since charge exchange and electron impact are not included as losses. 

490 In the second case (raii) , we derived a lower limit to the lifetimes from Table 2 of Fleshman 

491 et al. [2010b] by summing the additional losses due to charge exchange and electron 

492 impact, finding: H, O, OH, H2O = 0.4, 0.4, 0.2, 0.03 years, respectively. Notice in 

493 particular the drastically different times scales for H and O, where including the additional 

494 sinks reduce the size of the H cloud by a factor of 40/0.4 — 100, and the oxygen cloud 

495 by 14/0.4 = 35. This case represents an extreme limit, given that the neutrals spend 

496 almost all of their time orbiting outside of the Enceladus torus, where compared to photo- 

497 processes, the chances of charge exchange and electron impact are relatively unlikely. We 

498 mention, however, that Rymer et al. [2007, 2008] has shown that circulation patterns 

499 inside of 12 Rg at Saturn gives rise to 'butterfiy' hot electron pitch angle distributions, 

500 related to low temperature anisotropy {T±/T\\), on which proton field-aligned distributions 

501 depend [Sittler et al, 2008]. 
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502 The individual clouds (excluding hydrogen) were weighted by the stated time scales 

503 and totaled in Fig. 8b. When only losses to photodissociation/ionization are considered 
(''"phot)) the neutral cloud is dominated by oxygen, whose fate thus determines that of 

505 the neutral cloud. When charge exchange and electron impact are also included (rau), 

506 dissociated OH contributes significantly, driving the neutral cloud (ionized) percentage 

507 up, and the escape percentage down. We note that the neutral fates presented in Bagenal 

508 and Delamere [2011] (escape = 44%, ionized = 17%, absorbed = 39%) were based on an 

509 earlier version of our model which only included H2O. 

510 The particles that are neither absorbed nor lost by escape make up the neutral clouds. 

511 In the case where the cloud evolves for Tphot, oxygen and hydrogen dominate since they 

512 are far less hkely to be photoionized than are OH and H2O to be photodissociated. With 

513 charge exchange and electron impact included (ran), however, more oxygen and hydrogen 

514 are removed from the system, which then tends to favor a molecular OH-H2O cloud. In 

515 terms of total mass the same applies, although hydrogen accounts for only a few percent 

516 at most. We find that the total cloud mass is bounded between ^ 1 and 10 Mtons, for 

517 Tall and Tphot, respectively. 

518 It is worth pausing to re-emphasize that the system is in reality better represented by the 
Tphot case, from which all neutral clouds in this paper have been derived. The Taii case is 

520 strictly valid only for neutrals within the Enceladus torus, though reactions with electrons 

521 and protons may also prove important, as discussed above. What is illustrated, however, 

522 is that Saturn's magnetosphere is less oxygen-dominated than suggested by looking at 

523 losses from photo-processes alone. These results suggest that our oxygen abundances are 

524 somewhat overestimated, likely by less than a factor of two. 
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525 3.3.2. Neutral absorption 

526 The particles absorbed by Saturn and its rings are plotted by species and latitude in 

527 Fig. 9. In Fig. 9b, we see that most absorption comes from oxygen (74%), followed by H2O 

528 (11%), OH (9%), and finally by hydrogen (6%). Absorption is equally divided between 

529 Saturn and its rings except in the case of OH, where twice as much falls on Saturn's 

530 rings. This is because OH is largely produced by impact dissociation, which creates 

531 slower neutrals than does charge exchange, whereby in our model, H2O and oxygen arise 

532 exclusively. 

533 In Fig. 9b, absorption is plotted against Saturn's latitude. Because the model is sym- 

534 metric about the equator, the results apply to either hemisphere. Oxygen, water, and OH 

535 follow the same trends because they all originate from charge exchange (dissociated OH is 

536 slow and does not reach Saturn), and have been created from ions with similar velocity dis- 

537 tributions. Any second-order differences due to the velocity-dependence of the respective 

538 cross sections are not immediately apparent. Hydrogen, on the other hand, is produced 

539 entirely by dissociation in the model and exhibits a more uniform flux across Saturn. The 

540 explanation is that the velocity distribution from which hydrogen is produced is isotropic, 

541 whereas that which produces charge-exchanged neutrals is bi-Maxwellian (section 2.2.2). 

542 The fluxes shown in Fig. 9b are consistent with Hartogh et al. [2011], who modeled recent 

543 Herschel observations of Saturn's water torus and found an average flux of 6 x 10^ cm~^ 

544 s"^ for H2O -|- OH impinging on Saturn. 

4. Discussion 

545 Some useful conclusions can be drawn by further contrasting our results with Cassidy 

546 and Johnson [2010] (CIO). It is important that we flrst mention a profound difference 
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547 between our models. The model of CIO effectively carries out resonant charge exchange 

548 only, which does not chemically alter the neutral population; neutrals in their model are 

549 produced either directly from Enceladus or from subsequent dissociations. Neutrals in 

550 our model, on the other hand, originate from Enceladus (H2O). OH is then created via 

551 dissociation (as with CIO), but secondary O, OH, and H2O populations are created from 

552 H2O via charge exchange with the dense plume-fed Enceladus torus. The CIO model re- 

553 distributes neutrals around Saturn, while we redistribute and chemically re-assign neutral 

554 abundances by allowing for asymmetric charge exchanges. Thus, it may well be a coinci- 

555 dence that our models are similar in total abundance. While it may be difficult to compare 

556 our total abundances, the slope of our radial density profiles can be contrasted directly 

557 because our redistribution mechanisms (charge exchange and dissociation) are similar. 

558 Differences are due largely to ClO's inclusion of neutral collisions and our prescribing 

559 unique velocity-dependent charge exchange for each of the 0-, 0H-, and H20-producing 

560 reactions (reactions 3a-3c). 

561 Our neutral clouds are compared with CIO in Fig. 10. All of our clouds include contri- 

562 butions from charge exchange, but the H2O cloud is mostly comprised of water sourced 

563 directly from Enceladus (3.95 Rs), and OH includes the additional source from dissocia- 

564 tion. In the CIO model, the water molecules were spread due to neutral-neutral collisions, 

565 which explains our higher H2O densities near Enceladus's orbit (Fig. 10a). The slope of 

566 the oxygen profile agrees best with CIO because their charge exchange cross section most 

567 resembles our own (Eq. 16c). Our H2O profiles agree less, and our OH slopes, the least, 

568 due mainly to the strong effect that neutral collisions have on those more polar molecules. 

569 In particular, CIO used a much larger cross section for neutral collisions involving H2O 
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570 and OH [Teske et al, 2005] than for atomic oxygen [Bondi, 1964]. This helps to further 

571 explain our agreement with their oxygen profile since we exclude neutral-neutral collisions 

572 from our model altogether. We conclude that neutral-neutral collisions appear to play a 

573 less significant role with atomic species, such as oxygen and hydrogen. 

574 The column densities (Fig. 10b) are similar to CIO, who constrained their O and OH 

575 clouds with the most recent Cassini UVIS results of Melin et al. [2009]. Our oxygen 

576 density — as well as our total oxygen content (Fig. 10c) — is higher for two reasons. First, we 

577 use a larger cross section than does CIO for reaction 3c, and second, the clouds presented 

578 here have been limited only by photoionization. Charge exchange and electron impact are 

579 second order losses beyond 6 Rs, but including them would favorably reduce the oxygen 

580 content more than OH and H2O (section 3.3), bringing our models into better agreement. 

581 Our total H2O content is 4x less than CIO found (Fig. 10c). This is partly because we 

582 have subjected H2O molecules in the primary (plume-fed) neutral torus to the shortest 

583 lifetimes possible (section 2.1.1), whereas CIO tracks molecules that get kicked out of 

584 the densest plasma via neutral collisions, and thus survive longer, being less susceptible 

585 to both charge exchange and electron impact. That their total H2O content is higher 

586 than ours (Fig. 10c), does not contradict the fact that their H2O column density is lower; 

587 neutral-neutral collisions would spread out the torus, lowering the column density, while 

588 allowing neutrals to survive longer, increasing the total abundance. 

589 Our model would benefit by including the redistribution attributed to neutral collisions 

590 by allowing particles to interact in a direct simulation Monte Carlo (DSMC) model such 

591 as in CIO. Likewise, DSMC models would benefit by including charge exchange cross 
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592 sections specific to each reaction. Such models should also take into account asymmetric 

593 charge exchanges, which affects neutral cloud composition. 

594 The reactions modeled in this study were chosen in order to measure the effect of 

595 symmetric and asymmetric charge exchange at low velocities. Building upon our find- 

596 ings, future studies should include additional neutral-producing charge exchanges, such as 

597 0H+ + H2O ^ OH* + H2O+, H+ + H2O ^ H* + H2O+, and 0H+ + H2O ^ O* + H3O+, 

598 as well as dissociative recombination of H2O"'". 

5. Conclusions 

599 We have modeled low-velocity charge exchange from the point of view of the ions, 

600 allowing us to study the effects of velocity as well as gyrophase. With reactions 3a-3c, 

601 we have been able to offer an estimate on the size and shape of the neutral clouds at 

602 Saturn, while simultaneously exploring the sensitivity of the neutral clouds to a variety 

603 of velocity-dependent reactions. 

604 We have also re-visited the production of OH following H2O dissociation in the primary 

605 neutral torus. Previous models have used 1 km s~^ as the initial velocity for OH, while 

606 measurements suggest a range of speeds from 1 to 1.6 km s~^. In our model, the higher 

607 speed increases the range within which dissociation dominates neutral production from 9 

608 to 15 Rs. 

609 Additional findings are: 

610 (1.) Charge exchange cross sections that increase steeply at low speeds tend to produce 

611 neutral clouds more confined to the orbit of Enceladus, implying the most spreading for 

612 oxygen, moderate spreading for H2O, and the least for OH (Fig. 3). Accounting for 
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613 gyrophase doubles the local OH density within Encelacus's orbit, has ~ no effect on H2O, 

614 and decreases oxygen density by less than 10%. 

615 (2.) Enceladus is solely responsible for the creation of the neutral H2O torus via thermal 

616 ejection from its plumes. However, Saturn's neutral clouds are overwhelmingly produced 

617 by charge exchange and dissociation occurring throughout the torus (99%), and not near 

618 Enceladus itself (Fig. 5). 

619 (3.) We estimate that roughly half of all neutrals escape the system, with the remaining 

620 equally divided between absorption by the rings/planet and the neutral clouds (Fig. 8). 

621 Less than 50 kg s~^ is thus ionized and transported out of the system as plasma. This 

622 number is expected to represent an upper limit, given we have assumed that all particles 

623 forming the neutral clouds are ultimately ionized; a more accurate result would require 

624 modeling the detailed effects of charge exchange and neutral-neutral collisions within the 

625 neutral clouds. This estimate can be compared to Sittler et al. [2008] , whose Figs. 14 and 

626 17 give roughly {[N L'^]^J^+ / L'^) x mw+ /nransport ~ 3 X 10^^ x mw+/10^ s a; 10 kg s"^ at 

627 L = 10. 

628 (4.) Saturn's neutral cloud has a total mass of at least 1 Mton, but likely much closer 

629 to 10 Mtons. The primary plume-fed neutral torus (0.3 Mtons) is comprised entirely of 

630 water in our model, while the secondary neutral clouds are broken down into H (< 5%), 
O (< 82%), OH (> 13%), and H2O (~ 1%). Atomic oxygen dominates the composition 

632 both because of a high production rate from charge exchange as well as a long lifetime 

633 against photoionization. Charge exchange and reactions with electrons favorably remove 

634 hydrogen and oxygen, but are secondary loss mechanisms throughout the majority of the 

635 magnetosphere. 
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636 (5.) Our model predicts fluxes on Saturn from charge exchange of ~ 6 x 10^ cm~^ s~^ 

637 for both OH and H2O (consistent with Herschel observations by Hartogh et al. [2011]), 

638 and oxygen is about 5x higher. Absorption is divided equally between Saturn and its 

639 rings (Fig. 9a). 

640 (6.) Our total neutral abundances are similar to Cassidy and Johnson [2010] (CIO) for 

641 both OH and H2O, and 4x higher for oxygen (Fig. 10). Differences in the slopes of our 

642 equatorial density profiles are in part due to our not including neutral-neutral collisions, 

643 while this fact appears to have no effect on the oxygen profile. On the other hand, CIO 

644 did not include the effects on neutral chemistry following asymmetric charge exchanges, 

645 nor did they use velocity-dependent cross sections particular to each reaction. Herschel 

646 observations by Hartogh et al. [2011] confirm the importance of neutral-neutral collisions 

647 for H2O, but if oxygen is the dominant neutral species in Saturn's magnetosphere, as 

648 our model predicts, neutral-neutral collisions may play a smaller role in Saturn's neutral 

649 cloud than previously expected. 

650 Given the effect on both the size and shape of the neutral clouds, we suggest that 

651 future neutral cloud models include charge exchange cross sections unique to each reaction. 

652 Asymmetric charge exchange also has an important effect on neutral chemistry that should 

653 be implemented. Regarding the ions' gyrophase, Monte Carlo models can account for 

654 its effect by using phase-dependent probability distributions. Finally, the range of OH 

655 velocities studied here should be considered when modeling dissociation. 

656 Moving forward, we plan to implement these suggestions into the neutral cloud model 

657 of CIO and to couple that model with the chemistry model of Fleshman et al. [2010b]. 

658 Constrained by Cassini plasma observations, the chemistry model uses ClO's neutrals as 
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659 input, and provides ion temperatures and densities throughout the magnetosphere (< 20 

660 Rs), which CIO in turn uses to update neutral densities. An improved understanding of 

661 two issues is planned: (1) Where does plasma transport become important? (2) What is 

662 the role of hot electrons with regard to ion-neutral chemistry inside 20 Rs? 

Appendix A: Collision probability 

663 The average number of collisions occurring during a time interval At is given by A = 

664 i/At, where u — nneutraiscl'^reO'yrei IS the local coUision frequency, assumed to be constant 

665 during At. Statistics are applied to determine if and how many reactions occur during At. 

666 The Poisson distribution [Zwillinger and Company, 1996; Reif , 1965] gives the probability 

667 of suffering exactly n collisions for a given A: 

/(n;A) = — . (Al) 

668 Notice that Eq. Al peaks at n = A if one treats n as a continuous variable. Summing Eq. 

669 Al discretely from n — k to n — cg gives the probability of suffering at least k collisions 

670 during At, 

^.(A) = e-^$:-. (A2) 

n=k 

671 Because Eq. Al is normalized {e~^ Z^^o '^"Z^' — e'^e^ — 1), Eq. A2 can be conveniently 

672 rewritten as 

k—l 

Pk{X)-l-e-'^^. (A3) 

n=0 

673 A random number (0 < < 1) is compared to each Pf. at each timestep. The largest k 

674 for which > N determines how many fast neutrals (collisions), k, are produced during 

675 At. 
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676 In practice, it is only necessary to compare to the first few Pk wlien A <^ 1, made evident 

677 by tfie leading terms in Eq. A2 for A; + 1 and k: 

n+i^ /(fc + l;A) _ A^+V(fc + l)! _ A 

678 Multiple collisions are thus increasingly unlikely when A <^ 1. In such cases, comparison 

679 with Pi — 1 — ~ A = ulst is sufficient. 
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Figure 1. Sketch of gyrating ions in the neutral frame with guiding centers moving along 
a prescribed flow field, shown here near Enceladus for scale. Warm ions {v± > fflow) move on 
trajectories that coil around themselves and do not reach zero relative velocity with respect to the 
neutrals at any point. Cool ions {v± < tiflow) essentially trace their guiding centers with 'snake- 
like' trajectories, and also do not obtain zero relative velocity. Fresh pick-up ions {v± ~ ffiow) 
do, however, obtain zero relative velocity at the cusps of their cycloidal trajectories. Neutrals 
produced by charge exchange (whose trajectories are indicated by the red lines) tend to be created 
with velocities at which the respective reaction rates peak (Fig. 2). 
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Figure 2. (a.) Cross sections for the reactions listed in the legend. Data for reactions 3a and 
3b are from Lishawa et al. [1990], and data for reaction 3c is from Li et al. [1995]. ^^^^'■^poi^t'^*^^ 
^grmmetnc^ and (j^^y™™'=*"= qj-q hypothetical fits applying to the OH*-producing reaction, and are 
explored in Fig. 4. Ions oscillate between and 36 km s~^in the Enceladus torus, (b.) Collision 
frequency, na{v)v, for a given density of nnjO = 10^ cm~'^ plotted over the same energy range. 
The collision frequency increases with energy in the oxygen-forming reaction, while the water- 
forming reaction is independent of energy and the OH-forming reaction declines with energy. 
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Figure 3. Neutral clouds produced by the reactions shown in Fig. 2. (a.) Oxygen is the 
most abundant because the cross section section is 10 x higher than than with O and OH. The 
lifetime of oxygen against photoionization is also much longer than the lifetime for either OH or 
H2O against photodissociation. (b.) Same as above, but normalized to peak. Oxygen shows the 
most spreading because reactants are produced with higher velocities (Fig. 2b), which expands 
the cloud. The same trend holds with H2O and OH, where OH tends to be created with the 

1 ; 1 • /TT o extrapolated \ 

lowest velocities (rig. 2, ctoh )• 
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Figure 4. Neutral OH clouds produced from three hypothetical charge exchange cross sections: 

extrapolated symmetric j asymmetric /-i-i- / N extrapolated i j-i i • i j. j -j. 

ctqjj , ctqjj , and ctq^ (-rig- (,&•) ctoh produces the highest density 
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, the lowest) at Enceladus because of the creation of additional low-velocity particles. 



(b.) Same as above, but normalized to peak. The differences in density in the tail is not an 
indication of spreading, but rather further illustrates the deficiency in the peak density, going 



r extrapolated , asymmetric 

from ctqjj to ctqi^ 
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Plume versus global charge exchange 
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Figure 5. (a) Comparison between charge exchanged neutrals produced near the Enceladus 
plume and those produced from the neutral torus as a whole — in this case for oxygen, (b.) 
Though shown here for oxygen, all charge exchange reactions near the plume result in a cloud 
with less spreading than their global counterpart due to the imposed slowing of the plasma 
(and hence, the release of slower neutral products) near the plume in response to mass-loading 
[Fleshman et ai, 2010a]. 
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Figure 6. OH clouds produced from charge exchange and high- and low-speed dissociation. 
Dissociation dominates neutral cloud production inside 9-15 Rs, at which point charge exchange 
becomes the dominant contributor. 
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Figure 7. Neutral cloud densities in the r-z plane, (a.) Hydrogen produced purely from H2O 
dissociation, (b.) Oxygen produced purely from charge exchange (reaction 3c). (c.) Hydroxyl 
produced from the combination of charge exchange and dissociation. Dissociation dominates 
inward of 9-15 Rs along the equator, while charge exchange (reaction 3b) tenuously fills the 
magnetosphere elsewhere, (d.) Water produced entirely by charge exchange (reaction 3a). (e.) 
Dense torus fed directly by the Enceladus plumes (section 2.1.2). 
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Figure 8. The fates of neutrals in our model along with the results from other models, (a.) 
Dissociation produces low- velocity neutrals and OH is thus not likely to escape or to be absorbed. 
Conversely, dissociation also produces hydrogen which largely leaves the system, (b.) The results 
of J06 [Johnson et ai, 2006], J07 [Jurac and Richardson, 2007], and CIO [Cassidy and Johnson, 
2010], along with our own weighted totals (excluding hydrogen; see section 3.3.1). In the case 
of Tphot; the lifetime of the cloud is determined by photoionization/dissociation only, whereas 
with Tall, we limit the lifetimes by also including electron impact and charge exchange. These 
limiting cases bound the previous studies, except that CIO has more absorption attributed to 
neutral-neutral collisions. 
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Figure 9. (a.) Neutrals absorbed by Saturn, plotted by species. Partitions with horizontal 
lines indicate percentages absorbed by Saturn's rings, (b.) Neutral flux on Saturn as a function of 
latitude. Neutrals produced by charge exchange (H2O, OH, and O) peak in flux at low latitudes 
due to the nature of the ion distributions from which they orginate, which have initial velocity 
vectors predominantly in the ring plane. Conversely, hydrogen flux is constant across Saturn 
because it originates from dissociation, whose velocity distribution is prescribed as isotropic. 
Note that OH produced by dissociation is not energetic enough to reach Saturn. 
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Figure 10. (a.) Total neutral clouds from our model, compared with Cassidy and Johnson 
[2010] (CIO). All clouds include contributions from charge exchange (reactions 3a-3c), while H2O 
is largely comprised of water sourced directly from Enceladus, and OH includes contributions 
from dissociation. The cloud densities are limited by photodissociation for OH and H2O and 
by photoionization for O. Including charge exchange as a loss for cloud neutrals would reduce 
the lifetime for O more than for either the OH or H2O, and would lower the relative oxygen 
abundance accordingly, (b.) Equatorial column densities found by integrating the plotted equa- 
torial densities. The H2O column density is similar to CIO, despite their having a very different 
radial distributions, (c.) Total neutral cloud content. Our total H2O content is less than CIO 
found, while our H2O column density is higher because our H2O cloud is not subjected to neutral 
collisions and is thus more confined. 
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